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The perturbative calculation of the effect of fluctuations in the nonzero frequency modes of a 
weakly interacting Bose gas on the condensation temperature is reviewed. These dynamic modes, 
discarded in most of the recent studies, have a temperature-induced energy gap that allows for a 
perturbative approach. The simple, yet powerful algorithm used to calculate the effect in a high- 
temperature expansion in conjunction with zeta function regularization of infrared divergences is 
explained in detail. The algorithm is shown to be reliable by demonstrating that it reproduces 
known results for a series of examples. With two-loop contributions properly included, the dynamic 

■ modes are seen not to lead to a shift in the condensation temperature, thus revising our earlier 

' finding obtained at one loop. 

(N ■ 
G : 

a ■ 

Hi ' 

I I. INTRODUCTION 

(N ; 

The past five years or so have seen numerous studies aimed to determine the condensation temperature of a weakly 
interacting Bose gas (for references and critical discussions of the various approaches, see Ref. 0). These studies, 
driven by experiments on Bosc-Einstein condensation (BEC), were based on a variety of analytic, computational, and 
experimental methods. With only a few exceptions, however, the theoretical studies took the same starting point 
I of considering only the static part of the quantum field describing the bosons. The rationale is that the dynamic 
t-H ■ modes have a temperature-induced energy gap, making them irrelevant compared to the static mode, for which that 
gap vanishes, when studying the long-distance critical behavior of the system close to the phase transition. With the 
gapped thermal modes decoupled, the microscopic quantum theory reduces to a static, classical theory, assumed to 
properly describe this thermal phase transition in equilibrium. Universal properties, such as the critical exponents 
can indeed be computed from this naive theory, as is well known. 

Critical temperatures are, however, not universal. Although phase transitions occurring in different systems may 
display the same long-distance behavior characterized by the same values of the critical exponents, they in general 
*^ ■ take place at different temperatures. For example, BEC in a relativistic free Bose gas is in the same universality class 
as BEC in its nonrelativistic counterpart, yet the condensation temperatures differ in both systems Q- This difference 
arises because the critical temperatures of these noninteracting systems are determined entirely by the dynamic modes 
J-^ ' (see Sec. Ill B|l whose spectrum differ for the relativistic and nonrelativistic Bose gases. The universal properties, on 
^ , the other hand, are determined by the static modes, which arc similar for both gases. Whereas universal quantities 
are largely independent of the dynamics and details of the underlying theory, nonuniversal quantities in principle 
depend on them. The critical temperatures of lattice models are even sensitive to such short-distance details as the 
^ ■ precise lattice used. 

The proper way to proceed when calculating nonuniversal properties is well established. One has to go down one 
level and start from the microscopic model describing the quantum system under consideration and perturbatively 
integrate out the irrelevant degrees of freedom (the dynamic modes when considering BEC). One thus obtains an 
effective theory in terms of the remaining, relevant degree of freedom (the static mode in case of BEC), with the 
coefficients being determined by the irrelevant degrees of freedom that have been integrated out. This effective field 
program was carried out for the problem at hand almost a decade ago 3j , with the conclusion that the dynamic modes 
lead to a shift in the condensation temperature of a weakly interacting Bose gas. The purpose of this paper is to 
revise that finding. Namely 0, the shift obtained at one loop is canceled by a two-loop contribution not considered 
in Rcf. 3J. There, only the leading contribution in a high-temperature expansion of the two-loop Feynman diagrams 
was considered, whereas the canceling terms appear in the next order. To show that it is this omission and not, as 
has recently been suggested 0, the method used that led to the incorrect conclusion, we clarify and justify in detail 
the perturbative approach of Ref. 3]. To appreciate the approach, we first recall some essentials of finite-temperature 
field theory 
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A. Finite Temperature 



The properties of a quantum system at finite temperature can be studied |y] by going over to imaginary time 
t — > — it, with t restricted to the interval < t < hp, where (3 = 1/k^T is the inverse temperature. The time 
dimension thus becomes compactified and, consequently, the energy variable hkg discrete k$ — > icu n , with lo u the 
so-called Matsubara frequencies, 

. J 27r/3~ 1 n for bosons 

" = \ TT/3- 1 (2n + 1) for fcrmions, ( > 

where n is a (positive or negative) integer. The difference between the bosonic and fermionic frequencies arises 
because a bosonic field 0(r, x) satisfies periodic boundary conditions 4>(h(3, x) = 0(0, x), whereas a fermionic field 
satisfies antiperiodic ones ip(hf3,x.) = — ?/>(0,x). Being (anti)periodic, the fields can be expanded in a Fourier series as 
(with a similar equation for tp): 

n— — oo x ' 

where </>„(k) := 0(w„,k). 

An integral over the energy variable hko at the absolute zero of temperature becomes a sum over the Matsubara 
frequencies at finite temperature: 

^(M-^E^-")' (3) 

where g is an arbitrary function. 

The frequency sums encountered below are typically evaluated using the identity 0] 
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or because of antisymmetry 



/3 ^ ftV, + z 2 2 



i + »(*). (5) 
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where n(z) is the Bose-Einstein distribution function 

n {z) = ^- i =l-- 1 - + O m . (6) 

The first term at the right hand of the identity @ gives the result for zero temperature, while the second gives the 
thermal contribution. Such a splitting is convenient because the pure quantum T = part can often be combined 
with the bare (zero-loop) theory. In a renormalizable theory, for example, the zero-temperature part leads to a 
renormalization of the parameters of the bare theory. These contributions are then combined with the bare theory to 
dress its parameters. 

In this paper, the thermal contributions to the grand canonical thermodynamic potential f2 of various systems are 
expanded in a high-temperature series. According to standard statistical mechanics, f2 is related to the partition 
function Z via 

Q. = -/3- 1 hi(Z). (7) 

Physically, this potential determines the thermodynamic pressure p through Q = —pV, with V the volume of the 
system. High-tcmpcraturc expansions are an important computational tool for studying interacting theories where 
usually no exact results can be derived as for free theories. The standard approach is to use the identity backwards 
and to replace the Bose-Einstein distribution function again with an expression involving a frequency sum: 

t \ 1 12 ^ — ^ z , . 

n(z) = ^-2 + ^^^T^' (8) 

n—l n 
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where in the first term at the right hand, the n — contribution is isolated. The second term cancels the T = 
contribution contained in the sum. These two terms correspond to the first terms in a Laurent expansion of the 
Bose-Einstein distribution function. 

The n = contribution, stemming from the </>o(k) mode in the Fourier series @, is of special interest to us. This 
zero- frequency mode, specific to bosonic systems, is a thermal mode for which the temperature-induced energy gap 
2irf3~ 1 n in the denominator at the left hand of Eq. (@J vanishes. It is related to the static field <!> alluded to above as 
follows: 



where use is made of the Fourier series (J2J. This classical, zero- frequency mode is the order parameter of BEC in a 
weakly interacting Bose gas. It develops a nonzero value in the ordered state where the U(l) symmetry (j> — > exp(i#)0 
(6 being a constant transformation parameter) is spontaneously broken, with |<I>| 2 determining the average number 
density of condensed particles. That is, a nonzero value of the order parameter (a classical object) implies the 
presence of a condensate (a quantum phenomenon). Being static, this mode is unaware of the time dimension and 
only notes the space dimensions. Compared to the dynamic modes n (k), with n ^ 0, it lives in one dimension less 
and its contributions are accordingly different in form from those of the gapped thermal modes, being typical for a 
d-dimensional rather than a (d + l)-dimensional theory. 

The following hierarchy can be identified, with the pure quantum zero-temperature modes living in d+1 noncompact 
spacetime dimensions, the gapped thermal modes also living in d + 1 spacetime dimensions, but with a compactified 
time dimension, and finally the zero-frequency mode living in d space dimensions, which can be understood as 
(d + l)-dimcnsional spacetime with the compact time dimension shrank to a point by letting (3 — > 0. Below, the 
zero- frequency mode contributions are labeled by a superscript to indicate that they correspond to = 0, while 
the zero-temperature (/? = oo) contributions are labeled by a superscript oo. The contributions described by the 
remaining terms in the Laurent expansion |JS} , finally, are labeled by a superscript {3. We will, somewhat loosely, refer 
to these last contributions as due to the gapped thermal modes, although they also encode the T = contributions, 
that are explicitly subtracted by the second term in the Laurent expansion. 

Zero-frequency modes can produce nonanalytic behavior, depending on the number of space dimensions. Dolan and 
Jackiw in their study of thermal phase transitions in relativistic quantum systems took a pragmatic approach toward 
these nonanalytic terms induced by such modes and simply discarded them when estimating critical temperatures. 
In light of the discussion above, this approach is justified as the critical temperature is to first order determined by 
the vanishing of the quadratic term in the effective theory Featuring the zero- frequency mode $, that theory is 
obtained by integrating out the gapped thermal modes, which give only analytic contributions. It is only in the next 
stage, where fluctuations in $ are studied, that nonanalytic contributions can appear. 



In this paper, regularization techniques routinely applied to handle ultraviolet divergent integrals are used to stem 
infrared divergences instead. Usually, such divergences are regularized by introducing an infrared cutoff. This type 
of regularization forms the infrared counterpart of the momentum cutoff regularization of ultraviolet divergences, 
where integrals are made finite by integrating only up to a finite momentum. Besides using momentum cutoff 
regularization, ultraviolet divergences can also be rendered finite by using dimensional regularization, where the 
integrals are generalized to arbitrary dimension, or analytic regularization, where powers appearing in the integrands 
are generalized to take arbitrary values. In the following, these techniques rather than introducing an infrared cutoff 
are used to regularize infrared divergences. 

The infrared divergences in this paper arise when expanding thermodynamic potentials in a high-temperature series, 
using the algorithm of Ref . 0] . The algorithm differs from the standard high-temperature expansion and focuses on 
the analytic terms needed for the study of the critical properties of an interacting system. A series of examples is 
given to show in detail that this high-temperature expansion in conjunction with the regularization techniques used 
is reliable, by demonstrating that known results are recovered. The examples include a nonrclativistic (Sees. Ill Al 
and III B| ) and a relativistic fSec. Ill Ojl free Bose gas, as well as their fermionic counterparts fSecs. ITTl AHIII C|) . and a 
BCS superconductor fSec. IIIID|l . In Sec. IIV Al fluctuations in the nonzero frequency modes of a weakly interacting 
Bose gas are considered by integrating out these gapped thermal modes in a loop expansion. In particular their effect 
on the condensation temperature is studied using the algorithm advocated here. As mentioned at the beginning of 
the Introduction, although irrelevant for the critical behavior itself, these high-energy modes can in principle affect 
the value of the condensation temperature, which-unlike the critical exponents -is nonuniversal and depends on the 
details of the underlying quantum theory. Owing to the temperature-induced energy gaps, the contributions due to 




(9) 
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the nonzero frequency modes can be calculated in perturbation theory. It is shown that the shift in the condensation 
temperature found in Ref. Q at one loop is canceled by two-loop contributions Q , not only in d — 3 but in arbitrary 
2 < d < 4 (Sec. IIV Bjl . The advantage of considering arbitrary dimensions is that BEC in two dimensions can be 
studied by taking the limit d — > 2. A discussion of some recent studies on the subject is given in Sec. IIV CI followed 
by conclusions in the last section (Sec.EJ. 



II. BOSONIC INTEGRALS 



A. Bose-Einstein Condensation 



As first bosonic system, a free Bose gas above the condensation temperature is considered. Being noninteracting, 
the field 4> n (k) describing the bosons can be integrated out exactly, leading to the thermodynamic potential f2 Q 

n a ! ~ r d d k , / . t n 2 k 2 \ 

n— — oo \ / \ / 

apart from irrelevant constant terms independent of /3 and /i. Here, m and A are the mass and the de Broglie thermal 
wavelength A = y/ '2irh 2 f3 j 'm of the particles, and /i < is the chemical potential. The argument of the logarithm is 
recognized as the kernel of the time-dependent Schrodinger equation in the imaginary-time formalism, describing a 
particle of mass m at finite temperature in a constant background potential —\x. As indicated by the single set of 
frequency sum and momentum integrals, Eq. (|10|l is a one-loop result that, because the system is noninteracting, is 
exact. The momentum integral is best carried out by using the Schwinger propertime representation of the logarithm 
and Poisson's summation formula 

oo oo 

E c2nina = E s ( a - w ) ( n ) 

n — — oo w— — oo 

to replace the summation over n by one over w. The integrals over the Schwinger propertime parameter and the 
momentum variable are then easily carried out to yield the familiar fugacity series 

o i 

w — l 

where t = d/2 + 1 and a = — f3fi > 0. In Eq. I|12[l. the singular zero-temperature contribution corresponding to 
it) = is omitted. A more careful derivation of the thermodynamic potential IjlOf) reveals the presence of an extra 
convergence factor exp(iw„?7) typical for nonrelativistic theories, where 77 is a small positive parameter that is taken 
to zero only after the frequency sum is evaluated . This convergence factor suppresses the T = contribution (see 
below). Physically, the series (|12|l denotes the sum over closed worldlines that wind w times around the imaginary 
time axis jlOj . The first factor in the summand measures the configurational entropy of the loops, while the second, 
with a being proportional to the worldline tension, is a Boltzmann factor, weighing loops according to their length 
measured by w. Loops with w > 1 correspond to exchange rings of w particles which are cyclically permuted after 
an imaginary time h/3, as they appear in Feynman's theory of the lambda transition in superfluid 4 He |Tlj . 

The expression (jlUfl can alternatively be evaluated by keeping the momentum integral to the end and instead of the 
Poisson formula using the identity For nonrelativistic free theories, the summand there includes the convergence 
factor exp(iw„T7) that suppresses the T = part 0. The thermodynamic potential then takes the well-known form 
1 

y = -jgF*/2+i(<*), (13) 



where 



2 r°° n 2T - 1 

FAa) = W)l dq ^~i (14) 



with q the dimensionless loop variable defined by q 2 = f3h 2 k 2 /2m. 
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On comparing the two alternative evaluations of the thermodynamic potential it follows that 

Y w- T e- aw = -— / dq-i- -, (15) 

where in the first representation the momentum integral has been carried out, while in the second it still has to be 
performed. The equivalence of both results is easily established by expanding the integrand in a geometric series and 
using the integral representation of the Gamma function. 



B. High- Temperature Expansion 

Next, the function F T (a) is expanded in a high-temperature series. As mentioned in the Introduction, such an 
expansion is an important computational tool for interacting theories, where exact results are absent. The standard 
approach is to use the identity (@J backwards, i.e., the Laurent expansion © of the Bose-Einstein distribution function. 
The zero-frequency mode, corresponding to the first term in the Laurent expansion, gives the contribution 

o poo „2t-1 

F ^ = W)l d VT^ = QT ~ lr(1 - T) ' (16) 

which for d = 3 leads to a nonanalytic contribution of the form ~ c?l 2 . 

As an aside, this contribution due to the zero-frequency mode can also be obtained from the worldline loop repre- 
sentation by replacing the sum over w in F T (a) by an integral: 

°° />oo 

Vur T e- c ™ -> / dww- T c- aw = a T - 1 T(l-T) = F°(a). (17) 

w = l J ° 

The zero- frequency mode c/>o(k) is important for our purposes as it determines the critical properties of the free 
Bose gas. To illustrate this, we differentiate the zero-frequency mode contribution twice with respect to the chemical 
potential to obtain the known result that the compressibility k of an ideal Bose gas diverges ass~ (— /i) T ~ 3 when 
the condensation temperature T is approached from above, where /i vanishes as ji <~ — (T — T ) 2 ^ d ~ 2 ' ) . For d = 3, 
this implies that 

^(-^-^-(T-Tq)- 1 , (18) 

and gives as value for the specific- heat critical exponent a co = — 1 |l2j| . Describing a free theory, this one- loop result 
is exact. 

Despite being noninteracting, a free Bose gas is not in the universality class of the Gaussian model, but in that 
of the spherical model [12J. The nontrival exponents derive from the constraint that the total number of particles 
be fixed. Without this constraint, the chemical potential would vanish instead as fi ^ Tq — T irrespective of the 
dimension. This would result in Gaussian exponents describing BEC in a free Bose gas at constant pressure 0. In 
four space dimensions, corresponding to the upper critical dimension, and above (i.e., d > 4), the exponents of the 
spherical and Gaussian model coincide. 

As mentioned before, the zero-temperature contribution vanishes because of the presence of a convergence factor in 
nonrelativistic free theories. More in line with the present approach, it vanishes even without including this factor, 

1 f°° 

-F™(a) = — d^-^0, (19) 



r(r) 

because integrals over polynomials give zero within the dimensional regularization scheme [l3l | . 

The rest of the terms in the high-temperature expansion, corresponding to the last terms in the Laurent expansion 
arc analytic. These contributions are obtained in the standard approach by expanding the summand in a high- 
temperature Taylor series before carrying out the momentum integral and the frequency sum, with each term giving 
a Riemann zeta function 



oo 

CW =E^7- (20) 



m 

m— 1 
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Specifically, one finds Q 

= g^(-")'. <*> 

with a the dimcnsionlcss expansion parameter. The condensation temperature is determined by these contributions 
due to the gapped thermal modes. At that temperature, the chemical potential vanishes, so that the equation for the 
particle number density n, 

n = -vW (22) 

with F T (a) in Eq. (|13|l replaced with F^(a), gives the condensation temperature 

2nh 2 ( n \ 2/d . . 

In other words, although the critical exponents of BEC in a free Bose gas are determined entirely by the static, zero- 
frequency mode, the condensation temperature, being nonuniversal, is determined entirely by the dynamic, gapped 
thermal modes. 

It is natural to ask whether the series (|21(l can also be extracted directly from the integral representation given in 
Eq. I|14|) . where the frequency sum has already been carried out, without returning to a representation involving such 
a sum. To show that the answer is affirmative, we expand the integrand there in a high-temperature series. Each 
term thus generated contains an integral of the form 

— — I—— (24) 
q cpi- - 1 ; 

and derivatives thereof with respect to the parameter p (which is set to unity at the end). For t < 1, these loop integrals 
diverge in the infrared. We handle the divergences using zeta function regularization, by analytically continuing the 
following equation 

--£-^ = nmt)p-\ (25) 

q cPi — 1 

with t initially chosen large enough so that the integral converges, to arbitrary values of t. With this regularization 
scheme, the high-temperature expansion H21J) easily follows, 

2 ,00 2r-l - C(r _n 

W)l dq ^~i ^ ^ WTT) { - a) H[al (26) 

showing that the contributions due to the gapped thermal modes can indeed (when properly regularized) be extracted 
directly from the function F T (a) in Eq. (|14|l without the need to return to the representation involving a frequency 
sum. This result is quite general and will in the next subsection be illustrated for a relativistic system. 



C. Relativistic Free Bose Gas 



To show that the conclusions of the previous section by no means hinge on the form of the nonrelativistic spectrum, 
leading to the specific structure of the integrand in Eq. (fT4"j) . the same analysis for a relativistic free Bose gas is given 
in this section. The system is considered at zero chemical potential /i = 0. It has no phase transition then, but the 
algebra is simpler than with a chemical potential included [15j , without losing any of the generalities we wish to point 
out. 

The analog of expression UlUfl for the nonrelativistic Bose gas reads ^(| 

n— — oo v 7 

with c the speed of light. The argument of the logarithm is now recognized as the kernel of the Klein-Gordon equation 
in the imaginary-time formalism, describing a relativistic free scalar particle of mass m at finite temperature. Being 
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a noninteracting system, the one-loop result l|27|l is exact. With the help of the identity (J2J, the thermodynamic 
potential becomes 



n 
v 



d d k 



(27T) 

(hc)- d f3-( d+1 '> 

( 47r )d/2 



— \/ h 2 c 2 k 2 



/3 _1 ln f] 



r(d/2) 



r(2r; 



r(r) 



H T {f3mc 2 ) 



(28) 



where the new integration variable is q = /3Scfc and the thermal part of the thermodynamic potential is determined 
by the function 



H T {y) 



1 



dq 



-,2-r-l 



r(2r)7 ^Ty^eV^ -i' 



(29) 



with t = d/2 + 1. This function 0| is the relativistic analog of the function F T (a) in Eq. (|14|l . 

The high-temperature expansion of the function H T (y) in d = 3 was computed by Dolan and Jackiw Q- To 
deal with logarithmic divergences that arise because four spacetime dimensions corresponds to the upper critical 
dimension, dimensional rcgularization was used by considering d = 3 — e, with e taken to zero at the end. For 
convenience, r = 5/2 — e/2 was only used in the integrand, whereas e was set to zero from the start in the prefactor in 
Eq. l(2*9"|) . These authors used Eq. © backwards, i.e., Eq. (JHJ), expanded the summand in a Taylor series, integrated 
term by term, and finally carried out the remaining frequency sum in each term, with the result Q: 



H 5/2 (y) 



7T 

360 



90 



y 



48 



y 



l 

128 



In 



y 4 + o(y e ), 



(30) 



where 7 is Eulcr's constant and y — /3mc 2 the dimcnsionless expansion parameter. The nonanalytic cubic contribution 
is due to the zero-frequency mode represented by the first term in the Laurent expansion of the Bosc-Einstein 
distribution function in Eq. (|29Jl : 



HIM = 



1 



r(5) 



Aq^r- 



y- 



48 



y - 



(31) 



where analytic regularization is used to handle the ultraviolet divergence. The second term in the Laurent expansion 
of the Bosc-Einstein distribution function gives 



1 1 



dq- 



2 r(5) Jo x/^+y 



1 

128 



7 

12 



(32) 



while the gapped thermal modes yield 



Hl /2 {y) 



7T 

360 



96 



y 



1 

128 



ln(27r) + 7 



y 



0(y% 



(33) 



so that H 5 / 2 (y) = -H5/2 (y) ~ ^5/2 (y) + ^5/2 (y) ■ ^ n cx ^ ra minus sign is included in the left hand of Eq. 1321) , so that 
-ff°/ 2 is the T = contribution that has to be subtracted from H^(a) to arrive at the purely thermal contributions 
H r (a). 

We next rederive this result using the algorithm given above and extract it directly from the integral representation 
(|29|l . by expanding the integrand in a Taylor series. Applying zeta function rcgularization to handle the infrared 
divergences as before [see Eq. (I25f) ]. we obtain 



1 



dq 



Taylor 



TT 

360 



96 V 



768 y 



1 

128 



i + 2C'(0) + 7 -| 
e 3 



y 



24576 y 



C'(~6) 10 
1474560 y 



0(y 12 ), 



(34) 



where this series is carried out to higher order. To show that this expression coincides with H^ 2 (y), we differentiate 
the identity [13 



C(s) = 2 s tt s - 1 sin (±tts) r(l - s)C(l - s), 



(35) 



8 



to obtain for s a (positive or negative) even integer: 

C'(s) = 2 s - 1 ir s cos (±tts) T(l - s)C(l - s), (36) 

or specifically 

C'(0) = -iln(27r), C '(-2) = -^C(3), C(-4) = ~C(5), C(-6) = "^C(7). (37) 

It now immediately follows that Eq. (|34|) gives the contributions due to the gapped thermal modes, so that the results 
obtained in the standard manner and the one used here indeed coincide. 

Given that also the zero-frequency mode contribution as well as the T = part can be extracted from the integral 
representation l|29|) . by identifying them with the first two terms in the Laurent expansion of the Bose-Einstein 
distribution function, it follows that no backtracking to a representation involving a frequency sum is required, 
provided zeta function regularization is applied. Moreover, the simple, yet powerful algorithm discussed here is easily 
implemented using an algebraic manipulation program to generate large orders in a high-temperature expansion for 
a wide class of problems, including fcrmionic systems (with the obvious changes), as will now be demonstrated. 



III. FERMIONIC INTEGRALS 



A. Free Fermi Gas 



In this subsection, a nonrclativistic free, spinless Fermi gas is considered. The thermodynamic potential is given 
by the one-loop expression 11U|) for a Bose gas with an extra minus sign included to account for the antisymmetry 
of the fermionic field. Using Poisson's summation formula the thermodynamic potential is readily written as a 
fugacity series @ 



o i °° 

w—1 



where r = d/2 + 1 and a = —flu > again. As for a Bose gas, this series can be understood as representing a sum 
over closed worldlines winding around the imaginary time axis w times, with each winding accompanied by a factor 
— 1 typical for fermions. 

Alternatively, the thermodynamic potential can be evaluated using the identity [6j 

oo 

a t-— 1 —ihh) n + z 2 

n— — oo 

where ui n = ir/3~ 1 (2n + 1) are the fermionic Matsubara frequencies and f(z) is the Fermi-Dirac distribution function 

The first term at the right hand of the identity 1)39(1 gives the result for zero temperature, while the second gives the 
thermal contribution. Nonrclativistic free theories include again a convergence factor which suppresses the T = 
part. The thermodynamic potential then takes the same form (|13() as for a Bose gas with the function F T (a) replaced 
with the function 

GJa) = —— / dq—i- . (41) 

n ' r(r) Jo V 2 +« + i v 1 

The equivalence of the two representations, implying 

we = - 

w — 1 

is easily established 



V {-l) w+1 w~ T c~ aw = — - / dq^- , (42) 
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B. High- Temperature Expansion 

As for a bosonic system, the standard high-temperature series expansion of the function G T (a) is to use the identity 
(Kfflll backwards and to replace the Fermi-Dirac distribution function again with an expression involving a frequency 
sum: 

™ = l- It (43) 

where the first term at the right hand subtracts the T = part encoded in the sum. Unlike bosonic systems, fermionic 
systems do not have a zero-frequency mode as none of the Matsubara frequencies vanishes. 

Following the algorithm advocated here, we instead derive this expansion directly from the integral representation 
of the function G T (a) given in Eq. (|41|) without first returning to an expression involving a frequency sum. Expanding 
the integrand in Eq. (|41|) in a Taylor series, we encounter loop integrals of the form 

q cPi +1 

which for t < are infrared divergent. These divergences are handled again using zeta function regularization, where 
the following equation 



q ePi 



I =r(t)(l-2 1 -*)C(t)p- t , (45) 



with t initially chosen large enough so that the integral converges, is analytically continued to arbitrary values of t. 
The high-temperature expansion now follows as: 

2 / d Q Taylor y, 2l _ T+l) y = y (4g) 

r(r) Jo e? 2 +« + l ^ v >T(l + l) y ' tV ' K ' 



1=0 



This result can be checked by substituting Eq. (|42|l and differentiating the resulting equation arbitrary many times 
with respect to the dimcnsionlcss expansion parameter a and setting that parameter to zero at the end, leading to 

oo 

J2(-l) W+1 ™- t =0--2 1 - t )at), (47) 

w= 1 

which is an identity. 



C. Relativistic Free Fermi Gas 



As last noninteracting example, a relativistic free Fermi gas at zero chemical potential is considered, consisting 
of particles (but no antiparticles, which would lead to an extra degeneracy factor) with half-integer spin a . The 
thermodynamic potential is given by the expression l|27|l for the Bose gas with an extra spin multiplication factor 
and a minus sign, i.e., — (2a +1), included. With ui n now denoting fermionic Matsubara frequencies, the one-loop 
expression is readily cast in the form 



n 
v 



(2a + 1 



(hc)~ d /3- 



-(d+i) 



f 



r(d/2) 



(47r)<V 2 

where the thermal part is determined by the function 

1 f° 
lAy) = ^W)Jo 

For r = 5/2, it has the high-temperature expansion 

7tt 4 it 2 , 



dqq d - V? 2 + (Pmc 2 y 



T(2r) 

r(r) 



I T (ftmc 2 



dq 



J1t-\ 



\Jq 2 + y 2 c^ q2+y2 + 1 



2880 192 



y 



l 

128 



In 



(!) 



7"4 



y 4 + o(y e ), 



(48) 



(49) 



(50) 



10 



where because of the absence of a zero- frequency mode, all terms are analytic in the dimcnsionless expansion parameter 
y — f3mc 2 . 

Apart from a minus sign, the T = contribution is the same as that for a Bose gas given in Eq. (|32|) : 



1 1 r 
= 2f(5) Jo 



dq 



q 



\/q 2 + y 2 



l 

128 



1 



--M?l- 5 



7 



(51) 



where again an extra minus sign is included in the left hand of this definition, so that gives the T = contribution 
(and not its negative). 

Using the algorithm discussed here, we expand the integrand of Eq. I|49(l in a Taylor series to obtain as contribution 
due to the gapped thermal modes: 



dq 



,2t-1 



Taylor 7tT 4 



+ 1 



2880 192 



y 



128 



.I-2C / (0)- 7 -21n(2) + ^ 



768 v ; 24576 s v ' 



127 



1474560 



C(-6)y w + 0(y 12 ), (52) 



which, with the help of the identities lj37J) and together with TgJ 2 (y)> reproduces the result Ij50|l (expanded to higher 
order). 

These last two examples illustrate that the high-temperature expansion we use in conjunction with zeta function 
regularization of infrared divergences is also valid for fcrmionic systems. 



D. High- Temperature Expansion of BCS Theory 



In Ref. [3j, the BCS theory was considered at finite temperature using the approach advocated here. This example 
is included to point out a shift in approach when going from a noninteracting to an interacting system. The main 
difference with a noninteracting system is of course that one- loop results no longer are exact, meaning that not all 
the degrees of freedom can be integrated out in an exact manner to obtain the thermodynamic potential as a known 
function of the thermodynamic parameters a and (3. Instead, irrelevant degrees of freedom are integrated out (usually 
only perturbatively) to arrive at an effective theory of the Landau form expressed in terms of an order parameter 
describing the phase transition. The order parameter, whose value is zero in the symmetric, disordered state and 
nonzero in the ordered state with spontaneously broken symmetry, determines the critical properties of the system 
close to the phase transition, where it disappears. 

After linearizing the quartic interaction of the BCS theory by a suitable Hubbard-Stratonovich transformation, 
the fcrmionic degrees of freedom can be integrated out exactly in a one-loop calculation to yield for a uniform BCS 
superconductor: 



Oi(A) 



= -2 



d 3 fc 



V ~J (2tt) 3 

Here, the frequency sum has been carried out and 



(53) 



E(k) = y/?(k) + |A| 



(54) 



is the BCS spectrum with the order parameter A providing an energy gap at the Fermi surface. This spectrum 
replaces the excitation spectrum £(fc) = h 2 k 2 /2m — /i of the elementary fermionic excitations of mass m in the normal 
state, where the energy is measured relative to the chemical potential [i, which, contrary to a dilute free Fermi 
gas, is positive in the weak-coupling BCS limit and given by the Fermi energy ijl = h 2 kp/2m, with hkp the Fermi 
momentum. The prcfactor 2 in Eq. I|53|l arises because the fcrmions come in two species, with spin up and down. 
The true grand canonical potential depending only on the thermodynamic parameters a and /3, is obtained from 
f2(A) by integrating out A. This can only be done perturbatively, with the simplest approximation corresponding to 
the saddle point of the integral, where the order parameter is considered to be a nonfluctuating field. 
To the one-loop expression l|53() . the tree contribution is to be added 



n (A) 

V 



fJB 



(55) 
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where g-Q < is the (bare) coupling constant of the local BCS interaction term, representing the effective attraction 
between fermions. The pure quantum zero-temperature term in Eq. 1)53|) gives, among other contributions, a quadratic 
term which can be combined with the tree contributions to yield the renormalizcd coupling constant g, 



1 

fjB 



d 3 k 



1 



(2tt)3 i£(fc)r 

The integral diverges in the ultraviolet, to regularize it we introduce a momentum cutoff HA to obtain 

1 _ 1 m 
g 5b 2ir 2 h 2 



A + O(A ), 



(56) 



(57) 



where the irrelevant finite part of the integral is omitted. The weak-coupling BCS limit corresponds to taking g# — > 
from below. 

With the loop integration approximated by 



d 3 fc 
(2tt) 3 



(58) 



where f(0) = mfcp /2ir 2 h 2 is the density of states per spin degree of freedom at the Fermi level, the thermal part of the 
thermodynamic potential (|53[) becomes of a form encountered in a relativistic free Fermi gas , with the dimensionless 
integration variable q = /3£. The approximation consists of extending the range of the £ integration from —fx < £ < oo 
to — oo < £ < oo. The problem thus reduces to one in a single space dimension (d = 1), corresponding to the value 
r = 3/2. To handle logarithmic divergences that arise, we again dimensional regularize the integral and consider the 
problem in 1 — e dimensions instead. We proceed in the same way as before and rather than returning to an expression 
involving a frequency sum by using Eq. (|39() backwards, we expand the integrand in Eq. (|29() in a high-temperature 
series and apply zeta function regularization to handle infrared divergences. This gives as contributions due to the 
gapped thermal modes: 



d 3 fc 
(2^)3 



In (l + c-M 2 W+|AP) = /T MO) { y + 



2C'(0)+ 7 + 21n(2) 



y 2 + U'(-2)y i + 0(y 6 ) 



(59) 



Because the order parameter appears only in the dimensionless expansion parameter y = /3|A|, the high-temperature 
expansion is tantamount to a Landau expansion in |A|. Contact with the standard approach is established through 
the identities l|37(l and the substitution 



1 



In (flhiOT,) , 



(60) 



where the Debeye energy h/jj^), being a measure of the inverse lattice spacing, is the physical ultraviolet cutoff and the 
temperature is the relevant infrared scale. This correspondence between the pole 1/e of dimensional regularization 
and the logarithm ln(A) appearing in momentum cutoff regularization is commonly used in the context of quantum 
field theory, and can be verified explicitly here by working with the Debeye energy as a cutoff from the beginning 
instead of using dimensional regularization. Adding the tree and the one-loop contributions, we then find 



Oo(A) + gi(A) 
V 



In (Phuv) + 7-M*72) 



"(0)5 



|/3A| 2 -^C(3)|/3A| 4 



(61) 



featuring the renormalized coupling constant g. Being a static field, the order parameter is the zero-frequency mode of 
the problem at hand. The critical temperature T c is determined by the condition that the coefficient of the quadratic 
term changes sign, yielding the standard result pij 



knT r 



(62) 



with g < 0. Using this expression for the critical temperature, we can put Eq. (|61(l in the canonical form ^ 



no(A) + n l( A) = _ /rMo) r. 



(63) 



valid close to the critical temperature. 
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This illustrates first of all that our method of computing the contributions due to the gapped thermal modes in 
a high-temperature series expansion in conjunction with zeta function regularization of infrared divergences is also 
valid for an interacting system. In addition, it highlights the general approach to critical phenomena, which is to 
integrate out the irrelevant degrees of freedom (which include the gapped thermal and thus all fermionic modes) to 
determine the coefficients of the effective theory expressed in terms of the order parameter. Although irrelevant for 
the critical behavior itself, the gapped thermal modes can influence the value of the critical temperature through 
their contributions to the coefficient of the quadratic term. The critical temperature I)fi2(l . which is a one-loop result, 
is exclusively determined by the high-energy dynamic modes as none of the Matsubara frequencies of the fermionic 
fields that have been integrated out vanishes. Unlike critical exponents, the critical temperature is not universal and 
can depend on the details of these high-energy gapped modes. 



IV. CONDENSATION TEMPERATURE OF WEAKLY INTERACTING BOSE GAS 



A. Effective Theory 



We next consider the condensation temperature of a uniform, weakly interacting Bosc gas. Specifically, we are 
interested in how fluctuations in the irrelevant degrees of freedom influence this temperature. To determine this, we 
need, according to the effective field program (see the discussion in the preceding section), to calculate the coefficients 
of the effective theory expressed in terms of the order parameter by integrating out the irrelevant degrees of freedom. 
For BEC the irrelevant degrees of freedom are the gapped thermal modes, whose contributions arc conveniently 
computed using the algorithm advocated here. The order parameter $ of BEC in a weakly interacting Bosc gas, 
introduced in Eq. 0, denoting the expectation value of the quantum field 4> describing the bosons, is the zero- 
frequency mode. A nonzero $ implies the presence of a condensate, with |$| 2 giving the average number density 
"vac of condensed particles. This effective field prog ram is the /jeZrf-theoretic generalization of the effective classical 
potential introduced by Feynman and Kleinert [19j in their pathintegral study of thermal fluctuations via a final 
integral over the zero-frequency component of the path variable. (See the textbook |20j for extensive applications.) 
Here we are left with a final functional integral over the purely space-dependent field $(x). 

Whereas the condensation temperature of a free Bose gas is determined from within the normal, i.e., the high- 
temperature state, we determine T c for the weakly interacting gas instead from within the low-temperature state with 
a nonzero condensate |$|. 

The one-loop contribution to the thermodynamic potential of a weakly interacting uniform Bose gas is given by Q 



f d 3 k 



V J (27r) a 



c 



~PE(k)\ 



(64) 



where E(k) is the spectrum of the elementary excitation in the superfluid state expressed in terms of the order 
parameter 4> 



E(k) = \J (h 2 k 2 /2m - n B + 4 ffB |$| 2 ) 2 - 4g||$| 4 , (65) 

where and g& arc the bare parameters of the theory, with ge > the coupling constant of the repulsive two-particle 
contact interaction. The true thermodynamic potential Q(a, /3) is obtained by integrating out the remaining degree 
of freedom, which in a first approximation is done by evaluating the integral in the saddle point, corresponding to 
extremizing f2(<&). 

Written in the form (|65(l . the spectrum is the bosonic analog of the fermionic BCS spectrum l|54|) given in terms 
of the superconducting order parameter A. The spectrum l)65[l reduces to the T = Bogoliubov spectrum [2lJ when 
the lowest-order value |$| 2 (= n vac ) = /^b/2<7b for <& is inserted. This value is obtained from minimizing the tree 
contribution 

^ = - A ; B |<I>| 2 +.gB|$| 4 . (66) 

It is important to note that the chemical potential /is = 2gB"vac is positive here (so that the coefficient of the 
quadratic term is negative, leading to a potential of a form resembling a Mexican hat), setting a weakly interacting 
Bose gas apart from a free Bose gas which has a nonpositive chemical potential. 

The pure quantum zero-temperature term in Eq. (|64|) gives, among other contributions, two terms of the form 
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contained in the bare theory l|66[) . They can be combined to yield the renormalized parameters [2S 



H = mb - .9bA 3 



9 



9b 



2 r> 2 



(67) 
(68) 



with HA the momentum cutoff. In addition, it gives a nonanalytic contribution oc m 3 / 2 (c/|<I>| 2 ) 5 / 2 induced by the gapless 
Goldstone mode. This one-loop contribution, which is irrelevant for our present purposes and will subsequently 
be ignored, is typical for a gapless contribution in five spacetime dimensions. The effective number of spacctimc 
dimensions is five here because in nonrclativistic quantum theories, where time derivatives appear in combination 
with two space derivatives, the time dimension counts double compared to the (three) space dimensions. 

Since Eq. 164|) is already a one- loop result, it is consistent to this order to replace the bare parameters fiB,9B there 
with the (one-loop) renormalized ones. The renormalized coupling constant g is related to the s-wave scattering length 
a via IHHl 



2nh 2 a 



9 



(69) 



The zero-temperature relation between the renormalized, physical chemical potential and the average number density 
«vac = |3>| 2 of condensed particles then reads Q 



fi(T = 0) = 2gn v 



Antra 



(70) 



where in the last step we ignored the so-called depletion of the condensate [21j and replaced n va c with the total particle 
number density, which is justified because the term is already of order g. 

Next, the thermal part in Eq. I|64|) is expanded in a high-temperature series using the algorithm discussed above. 
Writing 



d 3 /c 



In 1 



32tt 1 

( 47r )3/2 ^3 J q 



dqq 2 In ( 1 — e 



-/3E(q) 



with q the same loop variable as used for a free Bose gas, defined by q — (3h k /2m, and 



0E(q) = yj(q* + a + Ayf 



Ay 2 



(71) 



(72) 



where y = (3g\<&\ 2 , we arrive at 



dqq 2 In ( 1 - e 



/3E(q)\ Taylor 



dq 



<f In 1 



4y) 



1 g V 

2 (ei 2 - 



1 



1 1 



2e^ 



-Ay 2 + O (f3 3 ) 



4 r (f ) C(f ) " C(|)(« + %) + |C(|)(« + 4y) 2 + + O (/3 3 ) 



(73) 



Contrary to the previous cases, we now have two dimcnsionless expansion parameters a and y, both being proportional 
to p. 

At two loop, we consider the class of bubble diagrams (see Fig. 4 of Ref. 0|) giving the following contributions to 
the thermodynamic potential due to the gapped thermal modes: 



9 



V 



(3/ 1 2 + 3/ 2 2 + 2/ 1 / 2 ), 



where I\ is the one-loop integral 



h 



d?k 1 h 2 k 2 /2m - n + 2g\<I>\ 



(2tt) 3 E(k) 



s 0E(k) _ 1 



(74) 



(75) 
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and a similar equation for I2 with 2g replaced with 6g in the numerator. Expanding the integrals in a high-temperature 
expansion, we obtain as two-loop contributions 

H^-l [C 2 (D -2C(i)C(|)(« + 4y)], (76) 

the leading term of which was considered in Ref. Q, but not the rest. 
Adding the contributions together, we arrive at 

Mtt =,-#fHi>r, (77) 

with 

co = ~ {C(|) - 2C 2 (§)<5 - C(|)« [1 ~ K(M + \C{\W} , (78) 
a ^-independent term, with expansion parameter 

5 = gp/\ 3 = a/A. (79) 

For <5 — > 0, Eq. (|78|l is analogous to the first contributions due to the gapped thermal modes of a free Bose gas given 
by Eq. JTHJ with Eq. (|2D>, albeit with a > now. Furthermore, 

-/3c 2 = [a + 4((M [l-4C(i)5], (80a) 
c 4 = g [1 - 12C(|)ff] • (80b) 

It is important to note that one cannot naively set <E> to zero in Eq. I|77|l [with VL2 ($) omitted], as was done by 
Baym et al. [| to argue that the zeta function regularization of infrared divergences used to obtain that result is 
flawed, and hope to retrieve the correct critical behavior of the compressibility, say, of a noninteracting Bose gas. 
As explained in Sec. B(l . universal properties are determined exclusively by the fluctuations in the zero- frequency 
mode. These fluctuations are, however, not considered in this section. We are interested here in the condensation 
temperature-a nonuniversal quantity. As argued in the Introduction, nonuniversal quantities in principle also depend 
on the gapped thermal modes. It is the effect of fluctuations in these dynamic modes that is summarized by Eq. J77J). 
More specifically, the coefficients given in Eqs. ifHIljl express the dressing of the parameters /1 and g by these thermal 
fluctuations. The one-loop corrections, first calculated in Ref. 0, have recently been rederived in Ref. p|, apparently 
in ignorance of the existing literature on the subject. 

The condensation temperature is determined by the vanishing of the coefficient C2, yielding 

a c = -4C(|)<5, (81) 

which gives T c in terms of /1. Experimentally, it is more realistic to consider the system at fixed particle number 
density and therefore take n rather than fi as independent variable. The particle number density is easily obtained 
from Eq. J77J at T = T c where $ vanishes: 



1 dfl 



L {C(|)[l-4C(|)5]-COc}. (82) 



a; 5 



T=T C 

Together with Eq. 1)81(1 this reduces to the free Bose gas expression: 

n = C(|)/A c 3 , (83) 

so that to this order the dynamic modes do not lead to a shift in the condensation temperature. This conclusion differs 
from our previous one in Ref. 0, where the two- loop contributions were not properly included 0. The additional 
two-loop contributions cancel the shift due to the gapped thermal modes found in Ref. 0] at one loop. 
With the expression l|83|) for n, the chemical potential at T c can be written as 

fx c = Agn = n, 84) 

771 

which is nothing but the Hugenholtz-Pines relation [2^| between the chemical potential and the self-energy at the 
condensation temperature where <f> disappears. It is identical to the zero-temperature relation l|70|) apart from a subtle 
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factor of 2 that can be understood by examining the Hugcnholtz-Pines relation, which is an exact result, at both 
temperatures. It follows from Eq. (|84J) that the chemical potential remains positive at the condensation temperature 
of a weakly interacting Bosc gas. This is in contrast to a free Bose gas (obtained by letting a — > 0), where the chemical 
potential vanishes when the condensation temperature is approached from above. Whereas the chemical potential in 
a free Bose gas remains zero all the way down to zero temperature, in a weakly interacting Bose gas it decreases from 
Agn at T c to 2gn at zero temperature. 

To justify the high-temperature expansion used, note that to lowest nontrivial order, the coefficient C2 is given by 
C2 = a + 4£(|)i5. The condition C2 = then gives the condensation temperature in terms of ji to lowest order as 



2wh 2 ( ii 



2/3 



fc B T = —^-3- , (85) 

which is a one-loop result. Since this temperature is large for g small, the high-temperature expansion is consistent 
with the weak-coupling assumption of perturbation theory. Using the Hugcnholtz-Pines relation (which is also satisfied 
at this lowest nontrivial order) to replace /j, with n, Eq. (|85|) takes the standard form (|23|) for a noninteracting Bosc 
gas, independent of g. As for a free Bose gas, this result is determined entirely by the nonzero frequency modes. 

Because a c cx 6 cx g, the coefficients (I80a() and (I80b|) of the effective theory are close to T c calculated up to the 
order g 2 . In the first coefficient (|78fl . the term cx 8 2 is missing. This term, given by the leading contribution in the 
high-temperature expansion of the third-order bubble diagram, is independent of the chemical potential and therefore 
irrelevant for our purposes. 

B. BEC in two dimensions 

For arbitrary dimension 2 < d < 4, Eqs. (|81() and l|83|) generalize to: 

a c = -4C(d/2)a, n = C(d/2)/A* (86) 

where now S = g(3/X d . The limiting case, d = 2, is special because ((d/2) diverges when d — > 2. To investigate this 
limit, we dimensional regularize the last equation by considering the problem in d = 2 + e, where 

2 1 

n =Af? (87) 

As in Eq. I|6(J|I . the pole in dimensional regularization can be connected with the logarithm appearing in the reg- 
ularization with a momentum cutoff, provided the ultraviolet momentum cutoff and the relevant infrared scale are 
identified. Here, they are given by the inverse range of the potential 1/a and the square root of the chemical potential, 
respectively |2(|, i.e., 

m ( -4=) = -M « k MMi/W)] (88) 



e \a^/mfj. J 2 \ln(l/87rna 2 )/ 2 

for ln[ln(l/87ma 2 )] >> 1. In deriving this use is made of the two-dimensional relation between the chemical potential 
and the particle number density at the condensation temperature, 

&Trh 2 n 
mm(l/8irna' ! ) 



With this correspondence, Eq. i|87(l leads to the well-known expression for the critical temperature |2fil l27j 

2nh 2 n , , 

k^T c = jrr. 90) 

to m|m(l/87mcr)J 

Together with the relation fj, c = Agn, Eq. (|89|) gives for the coupling constant 

2irh 2 
m\n(h 2 j mfia 2 ) 

To understand this result, recall that the time dimension counts double compared to the space dimensions in a 
nonrelativistic quantum theory. The quantum critical behavior of the theory was first investigated by Uzunov |28j |. 
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who showed that d = 2 (implying a total of four effective spacetime dimensions) corresponds to the upper critical 
dimension, below which a nontrivial infrared stable fixed point exists, with g > 0. For d > 2, the theory has a trivial 
infrared stable fixed point at g = 0. The logarithm in Eq. (|91J) arises because d = 2 separates these two different cases, 
with the chemical potential providing an infrared cutoff, so that a finite value for g and ultimately for T c obtains. In 
other words, the finite value for the condensation temperature in two dimensions hinges on the presence of quantum 
fluctuations. 

As an aside, we note that at zero temperature, the relation between the chemical potential and particle number 
density differs again a subtle factor of 2 from the relation ^ c = Agn we used to derive the critical temperature l|9U|l . 
As in d = 3, this factor is needed for the Hugenholtz-Pines relation to be satisfied at both temperatures. It implies 
that the an alog of Eq. (|89(l at zero temperature also differs a factor of 2, which is indeed what is found by Fisher and 
Hohenberg pf . 



C. Discussion 



Path-integral Monte Carlo simulations by Griiter, Ceperly, and Laloc |29f of the microscopic model describing a 
weakly interacting Bose gas restricted to relatively few (up to 216) particles in three dimensions found an increase in 
the condensation temperature linear in the scattering length, 

^-^ = Can^, (92) 
Jo 

with C = 0.34 ± 0.06. Since the dynamic modes do not contribute to the shift linear in a, the coefficient C can also 
be calculated in the static, classical theory. Recent Monte-Carlo simulations [3(j of that static theory gave larger 
values: C — 1.32 ± 0.02 and 1.29 ± 0.05, respectively. An intuitive understanding of the increase in T c was given in 
Ref. |29j |. where it was pointed out that a moderate repulsive interaction suppresses density fluctuations. This results 
in a more homogeneous system and facilitates the formation of large exchange rings necessary for BEC, which then 
takes place already at a higher temperature than in a free Bose gas. When the repulsive interaction increases further, 
the exchange is obstructed because it becomes more difficult for the particles to move. This leads to a lowering of 
the critical temperature as seen, for example, in liquid 4 He. A free gas with 4 He parameters at vapor pressure would 
have a condensation temperature of about 3 K, whereas liquid 4 He becomes superfluid at the lower temperature of 
2.17 K. 



; e sy i 

rings by slightly generalizing Eq. (|38|l for a free Bose gas |31j 



The full thermodynamic potential of an interacting Bose system can close to T c be written in terms of the exchange 

e gas |31j : 

oo 

J2^~ T e~ {a ~ ac)w , (93) 



n 

V 



oc 

w— 1 



where now the exponent r has the general form 



d 
D 



1, (94) 



with D the fractal dimension of the worldlines, which is D = 2 for a free Bose gas. The worldline tension determined 
by a — a c vanishes when the condensation temperature is reached from above as 

a-a c oc (T-Tc) 1 ^, (95) 

where a is a second exponent, which for a free Bose gas takes the value a = d/2 — 1 [sec above Eq. (|18(l ]. The two 
exponents a and r determine all the critical exponents characterizing the phase transition through scaling relations 
|31| . BEC in a weakly interacting Bose gas is in the same universality class as the lambda transition of liquid 4 He, so 
that the critical exponents describing BEC are accurately known. In going from a free Bose gas to strongly interacting 
4 He, the value of the fractal dimension D of the worldlines changes only slightly from D = 2 to D rj 1.96, with most 
of the interaction effects entering er, which changes from a = 1/2 to a sa 0.76 [21|. The value of the fractal dimension 
follows from the relation pjlj D — 2 — r), with ?] the Fisher exponent, determining the anomalous dimension of the 
order parameter 

For an ideal Bose gas, the value D = 2 derives from the quadratic form of the energy spectrum E(k) = h 2 k 2 /2m. 
A different value is obtained by modifying the free spectrum to E(k) ~ k D , with D ^ 2. It was noticed by Gunton 
and Buckingham 01 that an ideal gas with D = 3/2, so that E{k) ~ A; 3 / 2 , produces a value of the specific-heat 
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critical exponent a cc — close to the experimental value a cc m —0.01 for He |32j. However, this is a somewhat 
fortuitous and at the same time deceptive coincidence as other exponents come out incorrectly, and the actual value 
of the fractal dimension is close to 2, D « 1.96 j^- Describing the critical behavior of an interacting Bose gas, using 
quasiparticles with the spectrum 

E(k) ~ k 3 ' 2 , (96) 

as was done by Baym et al. can therefore not be justified. 

In addition to numerical studies of the static-classical theory, the coefficient C has also been estimated by analytical 
studies of that theory (for a summary, see Ref. such as the 1/N expansion |34j and variational perturbation theory 

M 

Although the dynamic modes do not, due to cancellations, contribute to the shift in the condensation temperature 
linear in a, they do determine the coefficient of the next term ~ a 2 n 2 / 3 lnfan 1 ' 3 ) in that shift 



V. CONCLUSIONS 



In this paper, the effect of fluctuations in the nonzero frequency modes of a weakly interacting Bose gas on the 
condensation temperature was studied in detail. The simple algorithm used to perturbatively calculate the effect in a 
high-temperature expansion in conjunction with zcta function rcgularization of infrared divergences was demonstrated 
to be reliable by showing that a host of known results arc recovered. The presence of temperature-induced energy 
gaps for these dynamic modes were argued to allow for a pcrturbative approach. It was shown that the shift in the 
condensation temperature of the form (|92J) with 

Cr = -|^«2.83, (97) 
we had obtained earlier at one loop 0] is canceled by two- loop contributions 0. 
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